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Abstract. Neutron stars and stellar-mass black holes are the remnants of massive stars, which 
ended their lives in supernova explosions. These exotic objects can only be studied in relatively 
rare cases. If they are interacting with close companions they become bright X-ray sources. If 
they are neutron stars, they may be detected as pulsars. Only a few hundred such systems are 
presently known in the Galaxy. However, there should be many more binaries with basically 
invisible compact objects in non-interacting binaries. 

Here we report the discovery of unseen compact companions to hot subdwarfs in close 
binary systems. Hot subdwarfs are evolved helium-core-burning stars that have lost most of 
their hydrogen envelopes, often due to binary interactions. Using high-resolution spectra and 
assuming tidal synchronisation of the subdwarfs, we were able to constrain the companion 
masses of 32 binaries. While most hot subdwarf binaries have white-dwarf or late-type main 
sequence companions, as predicted by binary evolution models, at least 5 % of the observed 
subdwarfs must have very massive companions: unusually heavy white dwarfs, neutron stars 
and, in some cases, even black holes. We present evolutionary models which show that such 
binaries can indeed form if the system has evolved through two common-envelope phases. This 
new connection between hot subdwarfs, which are numerous in the Galaxy, and massive compact 
objects may lead to a tremendous increase in the number of known neutron stars and black 
holes and shed some light on this dark population and its evolutionary link to the X-ray binary 
population. 



1. Introduction 

Neutron stars and stellar-mass black holes are the remnants of massive stars ending their lifes 
in supernova explosions. Detecting these exotic objects is possible when they are in a close 
orbit with another star. If matter is transferred from the companion star to the compact object, 
bright X-rays are emitted. Without ongoing mass transfer the companion remains invisible, 
but can be detected indirectly from the reflex motion of the visible star, which causes periodic 
variations of its radial velocity (RV). These variations are measureable via the Doppler effect 
from spectral line shifts. Stellar evolution models predict the existence of a hidden population 
of such compact objects. Subdwarf B or sdB stars are helium core burning stars of about half a 
solar mass with very thin hydrogen envelopes [Tj . A large fraction of sdB stars resides in close 
binaries [2], [3]. 



Because the components' separation in these systems is much less than the size of the subdwarf 
progenitor in its red-giant phase, these systems must have experienced a common-envelope and 
spiral- in phase [1], [5]. In such a scenario, two main-sequence stars of different masses evolve 
in a binary system. If the primary reaches the red-giant phase and fills its Roche lobe, mass 
is transferred to the companion star. When mass transfer is unstable, a common envelope is 
formed. Due to friction with the envelope, the two stellar cores spiral towards each other until 
enough orbital energy has been deposited within the envelope to eject it. The end product is a 
much closer system containing the core of the giant, which then may become an sdB star, and 
a main-sequence companion. When the latter reaches the red-giant branch, another common- 
envelope phase is possible and can lead to a close binary consisting of a white dwarf and an 
sdB star. In all known cases the companions are either white dwarfs or late-type main-sequence 
stars. 

2. Analysis method 

Since the spectra of close binary subdwarfs are mostly single-lined, they reveal no information 
about the orbital motion of the sdB stars' companions. However, in very close systems, the 
rotation of the sdB star is expected to be tidally locked to the orbital motion. This allows us 
to constrain the inclination angle of the system and to derive the companion mass (Figs. [Tj 
HJ). In order to derive the parameters, the mass of the sdB primary has to be known. Binary 
population synthesis models constrain the mass range for the sdB binaries in question. The 
mass distribution shows a sharp peak at about 0.47 Mq, which we adopted to derive the most 
probable companion masses. 

From the primaries orbital solution only the period P and the projected RV semi-amplitude 

K can be derived. The mass function f m = C ° m +M as)* = then provides a lower limit for 
the companion mass M comp for a given sdB mass M S( jb (Fig. 1). If the primary is synchronised 
(Fig. 2) the orbital period P equals the rotation period P ro t- 

• The stellar radius -R s dB is given by the mass radius relation R S( ib = \/M s ^G/g. 

• The surface gravity log g can be obtained by a quantitative spectral analysis (for details see 

El). 

• The rotational velocity v TOt = 2-7ri? sc iB jP can then be calculated and the projected rotational 
velocity v ro t sin i can be measured from the spectral line broadening of weak metal lines in 
the case of slow rotators [7] and of Balmer and helium lines for fast rotators [6] . 

Now the inclination angle i can easily be derived and the mass function can be solved for 
reasonable values of M s dB- 

The stars in our sample were observed with the high-resolution (R = 20 000 — 48 000) 
spectrographs UVES at the ESO VLT, HIRES at the Keck telescope, HRS at the Hobby Eberly 
Telescope (HET), FEROS at the ESO 2.2 m telescope and FOCES at the CAHA 1.5 m telescope. 

3. Nature of the unseen companions 

Today, 81 sdB stars in close binary systems have been studied and their orbital parameters have 
been derived (see the catalogue of Ritter & Kolb [8]). From 51 radial velocity variable sdBs in our 
sample, we selected 41 binaries with known orbital parameters or half of the known sample (see 
Fig. [3]). From these, 32 could be solved consistently under the assumption of synchronisation. 
In 9 cases the sdB primaries spin faster than synchronised. There are no spectral signatures of 
companions visible. Main sequence stars with masses higher than 0.45 Mq could therefore be 
excluded because of their high luminosities in comparison to the sdB stars. In this case spectral 
features of the cool secondary (e.g. Mg I lines at ~ 5 160 A) get visible in the spectra and a flux 
excess in the infrared appears, which can be detected using 2MASS photometry [9]. 



Figure 1. Schematic view of a single-lined binary system. 
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Figure 2. Schematic view of a single-lined binary system with synchronised rotation. 



Another possibility to detect M dwarf companions are reflection effects in the binary light 
curves. Some of our programme stars have already been checked for modulations in their light 
curves. Unfortunately, this method only works, if the binary inclination is high and the orbital 
period short enough. Synthetic light curve modelling of sdB+M systems shows that the expected 
amplitude of a reflection effect drops below lmmag for orbital periods longer than 0.5 d. This 
period therefore provides an upper limit for the detectability of reflection effects in sdB binaries 
from the ground. 

• In seven sdB binaries the companion has been identified as late M star. In addition to the 




Figure 3. Number of binaries plotted against the logarithm of their orbital periods. The solid 
blank histogram marks all known sdB binaries from the catalogue of Ritter & Kolb [8]. The 
shaded histogram marks the studied sample of 41 binaries with a focus on systems with rather 
short periods. 



derived companion masses, reflection effects are visible in their lightcurves. 

• Nine companions have to be white dwarfs, because of the derived masses and the absence of 
reflection in their lightcurves or flux contributions in their spectra. One sdB has a massive 
white dwarf companion and short orbital period. It qualifies as candidate for SNIa. 

• In seven cases this distinction cannot be made. The companions may be either late M stars 
or white dwarfs. 

• Eight systems have massive compact companions near or even exceeding the Chandrasekhar 
limit. Their companions therefore may be either unusually heavy white dwarfs, neutron 
stars or black holes (see Fig. 3). 



4. Orbital synchronisation of sdB binaries 

The assumption of orbital synchronisation in close binary sdBs is essential for our analysis. 
Theoretical models predict synchronisation to be established up to orbital periods of about half 
a day to two days [10], [TT]. But it has to be pointed out that tidal dissipation in radiative 
stellar envelopes is still poorly understood and the timescales of synchronisation differ by orders 
of magnitude depending on the models. Empirical evidence is therefore needed to constrain such 
models and to answer the question which close binary sdBs are synchronised. 

The timescale of the synchronisation process is highly dependent on the tidal force exerted 
by the companion. If the companion is very close and the orbital period therefore very short, 
synchronisation is established much faster than in binaries with longer orbital periods. If the 
sdB in a binary with given orbital period is proven to be synchronised, all other sdBs in binaries 
with shorter orbital periods should be synchronised as well. Although the timescales also scale 
with sdB radius and companion mass, the orbital period is dominating at first order. Most of 
the non-synchronised systems in our sample have orbital periods exceeding 1.3 d. 
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Figure 4. Mass distribution of the unseen companion stars. The companion mass is plotted 
against the total number of binaries under the assumption of canonical sdB mass 0.47 Mq. The 
solid vertical line marks the Chandrasekhar limit. The last bin at 7.5 Mq is a lower limit. A 
detail of the mass distribution is shown in the inlet. The shaded solid histogram shows the 
fraction of subdwarfs with confirmed compact companions, the shaded dashed one the detected 
M dwarf companions. The dashed vertical line marks the average WD mass. 



To study the influence of close companions on the rotational properties of sdBs we measured 
the projected rotational velocities of 49 single sdBs from the SPY survey [D]. We found that single 
sdBs are slow rotators spinning with an almost uniform rotational velocity of about 8kms _1 . A 
comparison with the v ro t sin i-distribution of our close binary sample revealed that the rotation of 
sdBs in close binary systems is clearly affected by the tidal influence of the invisible companions 
(see Fig. [5]). This alone does not prove that all sdBs in these binaries are synchronised, but it 
shows that the tidal influence of the companions is strong enough to significantly change the 
rotational properties of the sdBs. Although these stars may be not yet be synchronised, the 
synchronisation process is already at work. 

The parameters of the eclipsing binaries PG 1336-018 [H], HS 0705+6700 [13] and HWVir 
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Figure 5. Left panel The measured v ro tsmi of 49 single sdBs is plotted against relative 
fraction of stars as shaded histogram. The size fo the bins is given by the average error of the 
measurements. The blank histogram marks the expected uniform distribution of t7 ro tsini under 
the assumption of randomly oriented polar axes and the same rotational velocity v TOt = 8.3 kms -1 
for all stars. The dashed histogram shows this distribution for w rot = 9.0 kms -1 . The solid 
vertical line at f ro tsinz ~ 4.5 kms -1 marks the detection limit. All sdBs with lower v rot sini 
are stacked into the first bin (dotted histogram). Right panel The measured v r otsini of 51 RV 
variable sdBs is plotted in the same way. The size of the bins is given by the average error of 
the measurements and is therefore slightly different than in the left panel. The blank histogram 
marks the uniform distribution for w rot = 8.3 kms -1 . All sdBs with w rot sin i higher than 24 kms -1 
are stacked into the last bin. 



|14j derived from light curve analyses are consistent with the parameters derived with our 
method assuming synchronised orbits. This essentially means that the calculated u ro tsinz for 
synchronous rotation is consistent with the measured value. In eclipsing systems, all these 
parameters can be measured. This provides clear empirical evidence that at least the upper 
layers of the stellar envelope are synchronised to the orbital motion of the eclipsing sdB binaries 
in our sample. We therefore conclude that all sdBs in close binaries with orbital periods up to 
0.12 d should be synchronised as well. 

Two known sdBs clearly show ellipsoidal variations in their light curves with periods exactly 
half the orbital periods (KPD 0422+5421 [Td]; KPD 1930+2752 [IB], [§]). This alone is only 
an indication for tidal synchronisation. To really prove synchronisation it is necessary that the 
stellar parameters determined independently from the light curve analysis are consistent with 
a synchronised orbit. This is the case for KPD 0422+5421 as well as KPD 1930+2752. Both 
ellipsoidal variable systems have very short periods of ~ 0.1 d and high inclinations. Otherwise 
ellipsoidal variations are very hard to detect. 

This was possible in the case of the sdB+WD binary PG 0101+039 using a high precision 
light curve obtained with the MOST satellite [7]. We found a strong indication that the surface 
rotation of the sdB star PG 0101+039 is tidally locked to its orbit. Hence, other sdB stars in close 
binaries should also be synchronised if their orbital period is less than that of PG 0101+039 
(P = 0.567d). We conclude that tidally locked surface rotation is at least established in sdB 



binaries with orbital periods of less than half a day. 

An independent method to proof orbital synchronisation is provided by asteroseismology. Van 
Grootel et al. were able to reproduce the main pulsation modes of the short period pulsating sdB 
Feige 48 (P ~ 0.38 d), derived the surface rotation from the splitting of the modes and concluded 
that the subdwarf rotates synchronously [17] . Charpinet et al. reach a similar conclusion for the 
short period eclipsing binary PG 1336— 018 (P ~ 0.10 d) [18J. Asteroseismic analyses revealed 
that sdB binaries up to orbital periods of about 0.4 d are synchronised. We therefore conclude 
that all sdBs in close binaries with shorter periods should be synchronised as well. 

From the eight candidate sdB+NS/BH binaries found in our sample four have orbital periods 
shorter than 0.4 d where synchronisation could be proven by asteroseismology. Another two 
binaries have periods shorter than 0.6 d, the lower limit for synchronised rotation derived by 
binary light curve analysis. And the longest binary period is only 0.8 d. We therefore conclude 
that all these binaries should be synchronised and the derived parameters are therefore correct. 

5. Formation of sdB+NS/BH binaries 

The evolution that leads to such systems requires an initial binary, consisting of a primary star 
massive enough to produce a neutron star or black hole, and a companion of typically several 
solar masses. These systems experience two mass-transfer phases and one supernova explosion 
(see Fig. 4). The second mass-transfer phase had to be unstable, leading to a common-envelope 
and spiral-in phase. Population synthesis estimates suggest that of order 1 % of all hot subdwarfs 
should have neutron-star or black- hole companions [20], [19] . 

6. The fraction of NS/BH companions and selection effects 

This high fraction of massive compact companions (20 %) is in part caused by an interplay of 
several different selection effects (e.g. RV variability selection, bias on short period binaries, see 
Fig. E]). Especially the fact that most of the binaries with the highest derived companion masses 
have low inclination seems suspicious. It can be explained by the limitations of our analysis 
method. Nevertheless, high inclination systems should be more numerous than binaries seen at 
low inclination. The fact that no such system was found up to now (81 binaries with solved 
orbits are known), puts an upper limit of a few percent to the number fraction of sdB+NS/BH 
binaries. 

Both observational constraints and binary population synthesis models suggest a fraction of 
a few percent for sdBs with NS/BH companions. The fraction we found in our sample is too 
high to be consistent with these results even if selection effects are taken into account. It has 
to be pointed out that all of our candidate binaries have short orbital periods (0.2 — 0.8 d) and 
should therefore be most likely synchronised. Adopting the lowest possible mass for the sdBs 
(0.3 Mq , [1] , [5] ) only two of our candidate binaries (~ 5 %) have companion masses exceeding the 
Chandrasekhar limit. This number is consistent with theoretical and observational constraints. 
We therefore conclude that most sdBs in our candidate systems should have masses lower than 
canonical. 

7. HYPERMUCHFUSS survey 

The results presented here provide strong indications for the existence of sdB+NS/BH binaries, 
but rely on the assumption of orbital synchronisation and are seriously affected by selection 
effects. Most of the binaries with massive unseen companion are seen at low inclination. High 
inclination systems must exist as well and should be even more numerous. In this case a 
determination of the orbital parameters would be sufficient to prove their existence, because in 
binaries with such high RV amplitudes the lower limit for the companion mass (see Fig. 1) should 
already exceed the Chandrasekhar mass. The HYPERMUCHFUSS-survey has been launched 



Formation of subdwarfs 
with neutron stars 

Initial binary: Af, = 15iW Q , 
M 2 = 3M , P orb ~4yr. 




Formation of subdwarfs 
with black holes 

Initial binary: M, = 30 M , 
M 2 =5M , P orb ~10yr. 



First unstable mass transfer 
(when the primary is a red 
supergiant; case C transfer) 
leading to a common-envelope 
and spiral-in phase. 




First unstable mass transfer 
(when the primary is a red 
supergiant; case C transfer) 
leading to a common-envelope 
and spiral-in phase. 



Helium-star binary with 
M Be =5M Q and M 2 = 3 M Q , 
Porb =10d. 

Supernova (Type Ib/c) of the 
primary, producing a neutron 
star with M NS — 1.4 M Q and 
M2 = 3 M & in a highly eccentric 
orbit with P or b ~ 0.5 yr. 



Second unstable mass transfer 
(when the secondary is a red 
giant; case B mass transfer) after 
a brief X-ray binary phase. 




Helium-star binary with 

Miie = 10 M Q and M 2 = 5 M e , 

Porb = 20 d. 

Supernova (Type Ib/c) of the 
primary, producing a black hole 
with M B h ~ 4 M Q and 
Mi = 5 M @ in a highly eccentric 
orbit with P or b ~ 1 yr. 



Second unstable mass transfer 
(when the secondary is a red 
giant; case B mass transfer) after 
a brief X-ray binary phase. 
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P orb ~ 10 hr. 
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Helium-core burning hot 
subdwarf with 

Afadt, ~ 0.5 - 1 M Pj and a black 
hole in an ultracompact orbit 
with P orb ~ 10 hr. 



Figure 6. Schematic diagram of formation scenarios leading to hot subdwarf binaries with 
neutron-star (right hand panel) or black- hole (left hand panel) companions. 



in 2008 to search for sdBs with high RV variability and massive compact companions as well as 

the recently discovered class of hypervelocity stars (see Tillich et al. these proceedings). 
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